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y ! Abstract 

^ ' In constrast to discretized space-time approximations to contin- 

uum quantum field theories, discretized velocity space approxima- 
tions to continuum quantum field theories are investigated. A four- 
' momentum operator is given in terms of bare fermion-antifermion- 

boson creation and annihilation operators with discrete indices. In 
, continuum quantum field theories the fermion-antifermion creation 

and annihilation operators appear as bilinears in the four-momentum 
■ operator and generate a unitary algebra. When the number of modes 

op . range over only a finite number of values, the algebra is that associ- 

ated with the Lie algebra of U(2N). By keeping N finite (but arbitrary) 
problems due to an infinite Lorentz volume and to the creation of infi- 
^ . nite numbers of bare fermion-antifermion pairs are avoided. But even 

^ I with a finite number of modes, it is still possible to create an infinite 

number of bare bosons. We show how the full boson algebra arises as 
the contraction limit of another unitary algebra that restricts the num- 
ber of bare bosons in any mode to be finite. Generic properties of finite 
mode Hamiltonians are investigated, as are several simple models to 
see the rate of convergence of the boson contraction; the possibility of 
fine tuning the bare strong coupling constant is also briefly discussed. 
Classification: 81R10, 81R15, 81T05, 81T25, 81T27. 
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1 Introduction 



A central problem in liadronic physics is constructing nonperturbative so- 
lutions to QCD. One way to obtain nonperturbative solutions is to make 
space-time discrete and look for solutions on a space-time lattice. This is the 
lattice QCD program [1]. 

Another possibility is to discretize in momentum space; various groups 
have attempted to find Hamiltonian QCD solutions using instant and front 
form versions of discretized quantum field theory |2j. In this series of papers 
I explore quantum field theory on velocity grids, in the context of point 
form quantum field theory [5]. The motivation for this work is not only to 
obtain approximate nonperturbative solutions to Hamiltonian quantum field 
theories and in particular Hamiltonian QCD, but more generally to study the 
generic eigenvalue structure of trilinear and boson self-coupling interactions. 

In the point form [1] all interactions are in the four-momentum operator 
and Lorentz transformations are kinematic. Interactions are introduced via 
vertices, products of local free fields, which are integrated over the forward 
hyperboloid to give the interacting four-momentum operator. The natural 
variable that arises in point form is the four-velocity, the four-momentum 
divided by the bare mass of underlying constituents, and it is the space 
component of the four-velocity that is made discrete and finite. 

The four-momentum operator will be written as the sum of free and 
interacting four-momentum operators, = P^(/r) -|- P^{I). To guarantee 
the relativistic covariance of the theory, it is required that 



where [/a is the unitary operator representing the Lorentz transformation 
A. These "point form" equations [5], in which all of the interactions are in 
the four-momentum operator and the Lorentz transformations are kinematic, 
lead to the eigenvalue problem 



where p'^ is the four-momentum eigenvalue and |\I'p > the eigenvector of 
the four-momentum operator, which acts in generalized fermion-antifermion- 
boson Fock spaces. Then the physical vacuum and physical bound and scat- 
tering states should all arise as the appropriate solutions of the eigenvalue 








(1) 
(2) 



(3) 
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Eq.(3). What is unusual in Eq.(3) is that the momentum operator has inter- 
action terms. But since the momentum and energy operators commute and 
can be simultaneously diagonalized, they have common eigenvectors. One 
of the important properties of the point form is that the Lorentz generators 
have no interactions, so that global Lorentz transformations on operators 
and states are simple and explicit. 

Excluding boson self-interactions, all of the fundamental particle inter- 
actions have the form of bilinears in fermion and antifermion creation and 
annihilation operators times terms linear in boson creation and annihilation 
operators. For example QED is a theory bilinear in electron and positron 
creation and annihilation operators and linear in photon creation and anni- 
hilation operators. The well-known nucleon-antinucleon-meson interactions 
are of this form as are the weak interactions. These interactions differ of 
course in the way the fermions are coupled to the bosons, including the way 
in which internal symmetries are incorporated. For QCD, because of the 
SU{3)coior symmetry which generates gluon self-coupling terms, the gluon 
sector is no longer linear in creation and annihilation operators; this is also 
the case for the weak interactions and with gravitons; since gravitons carry 
energy and momentum, they also can couple to themselves. But even with 
bosonic self-interactions, the coupling of bosons to fermions is trilinear. 

If a^,b^ and denote respectively, fermion, antifermion and boson cre- 
ation operators, the aformentioned trilinear interactions can all be written 
as {a' + 6)(a + b^){c'' + c), while the "relativistic energy" terms are of the 
form a^a — 66^ + c^c. Written in this way the fermion-antifermion bilinears 
a^a, 66^, a^b^, and ba close to form a Lie algebra which is related to the Lie al- 
gebra of the unitary groups. Similarly the boson operators c^c, and c close 
to form a Lie algebra related to the semidirect product of unitary groups 
with the Heisenberg group. Then the aforementioned interactions can all be 
viewed as arising from an algebra of operators generated from these two Lie 
algebras. 

One of the main problems that arises in solving continuum field theory 
eigenvalue equations such as Eq.(3) is that the interacting four-momentum 
operator takes elements out of the generalized fermion-antifermion-boson 
Fock space. Difficulties arise in three ways, from the infinite Lorentz volume, 
from the possibility of creating infinite numbers of bare fermion-antifermion 
pairs, and from the possibility of creating infinite numbers of bare bosons. 
If the number of fermion-antifermion velocity modes is made finite, the first 
two kinds of problems can be avoided. As will be shown, for a finite number 
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of modes, N, there is an underlying fermion-antifermion symmetry generated 
by the Lie algebra of the group U(2N). With such a group structure it is 
possible to define an inductive limit as N goes to infinity [6] ; inductive limits 
will be studied in a following paper. Appendix B shows that for finite ap- 
proximations to trilinear interactions, the ground state energy varies linearly 
with the bare coupling constant (for large values of the coupling constant), 
with the slope being negative. 

But even with a finite number N of modes, it is still possible to have 
indefinitely large numbers of bare bosons in each mode. One of the main goals 
of this paper is to show that the algebra of bosonic operators, generated by 
the semidirect product of unitary with Heisenberg algebras, can be given as 
a contraction limit of another unitary algebra; before the contraction limit is 
taken, the number of bare bosons in any mode is finite. Then the eigenvalue 
problem, Eq.(3) becomes a problem in diagonalizing matrices. The maximum 
number of bare bosons is controlled by a number M, which, going to infinity 
as the contraction parameter goes to zero, gives the full boson algebra. 

For theories like Hamiltonian QCD there are also boson self-coupling 
terms. In that case the c^c terms in the free four-momentum operator are 
supplemented by self-energy terms. But the boson contraction limit is still 
valid with these self-coupling terms, only now the matrices to be diagonalized 
are more complicated. A simple one mode example is given in section 5 to see 
how fast the boson contraction converges; boson self-coupling terms are then 
added and the possibility of fine tuning the bare strong coupling constant is 
discussed. This section also analyzes a simple few mode isospin model. To 
further examine the nature of bosonic contraction. Appendix A discusses a 
simple exactly solvable bosonic model. 

Section 2 provides an introduction to point form relativity; starting with 
trilinear interactions, the existence of a unitary symmetry arising from fermion 
bilinears is discussed in section 3. Section 4 deals with bosons and boson con- 
tractions. 

2 Point Form Vertex Interactions 

The starting point in this paper for approximating the quantum field theory 
eigenvalue equation are the massive one-particle representations [7] of the 
Poincare group, with four-velocity states \v, a >, where v-v = 1, and the four- 
momentum is p = mv; a is the spin projection ranging between —j<a< j. 
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m and j are the mass and spin respectively. (For massless particles, v -v = 
and is treated separately) 

Under a Poincare transformation, 

Ua\v,cr> = e'^'°-\v,a> 

UA\v,a> = Y.l^^^'^' > K' jRw{v,A)), 

where a is a four-translation, D^, ^() is an SU{2) matrix element for spin j 
and Rw{v,A) is a Wigner rotation, an element of the rotation group S0{3) 
defined by 

Rw{v,A): = B-\Av)AB{v). 

B{v) is a boost, a Lorentz transformation satisfying v = B{v)v^'^^^, with 
^rest _ (1^0,0,0). Boosts are always canonical spin boosts [8J. 

From such one-particle representations, a many-body theory can be for- 
mulated by introducing creation and annihilation operators with the same 
transformation properties as one-particle states; these generate multiparticle 
states from the Fock vacuum: 

\v,a> = a'{v,a)\0>, (4) 
[a{v,a),a\v' ,a')]± = 2voS^{v - v')6^^^' (5) 
Uaaiv,a)U~' = e-'P-''aiv,a) (6) 

where the ± denotes commutator or anticommutator for bosons or fermions 
respectively. By using velocity rather than momentum variables, all opera- 
tors are dimensionless. 

Let ai,bi, Ck denote respectively bare fermion, antifermion, and boson an- 
nihilation operators where the indices stand for both space-time and internal 
variables such as charge or isospin. Then the free four-momentum operator 
can be written as 

P'^(/r) : = dvv^{ala, + bib, + Kclck), (8) 

where dv := — is the Lorentz invariant measure, k is a dimensionless relative 
bare boson mass parameter and m is a constant with the dimensions of 
mass; its value is determined by relating a physical mass such as the nucleon 
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mass to the dimensionless eigenvalue of the corresponding stable particle. 
The indices i and k contain both the continuous velocity modes and finite 
spin projection and internal symmetry degrees of freedom. Because of the 
transformation properties of the creation and annihilation operators inherited 
from the one particle states, the free four-momentum operator, as defined in 
Eq.(8), satisfies the point form equations (1) and (2) [5]. 

Interacting four-momentum operators are generated from vertices, prod- 
ucts of free field operators themselves made from creation and annihilation 
operators. Denote a vertex operator by V{x), where x is a space-time point. 
Then the interacting four-momentum operator is obtained by integrating the 
vertex operator over the forward hyperboloid: 



P^'(/) ■■ = g j (rx5{x ■ X - T')e{xo)x^'V{x) (9) 
= g dxx^V{x), 



with g a coupling constant. 

The vertex operator is required to be a scalar density under Poincare 
transformations and have locality properties. That is, 

UaV{x)U-' = V{x + a) (10) 
UAVix)Ux' = ViAx); (11) 
[Vix),Viy)] = (12) 

if {x — yY is spacelike. Here Ua = e"*^^-^''-''" is the free four-translation 
operator. 

Making use of the fact that if x and y are two time-like four-vectors 

with the same length, so that their difference is space-like, it follows that 

[P^, P'^] = 0. Also, from the Lorentz transformation properties of the vertex 

given in Eq.(ll), it follows that the interacting four-momentum operator 

transforms as a four-vector. Thus the interacting four-momentum operator 

also satisfies the point form equations (1) and (2). 

If the free four-translations are made infinitesimal, then [P'^{fr), -P^(/)] 
a - 



J dxx^ -!r-V {x) , SO that 



[P^{fr),P^^{I)]-[P^^{fr),P\I)] = g j dx{x^^~x^^)V{x) 



dxy dx 



= 0; 

[P^{fr) + P>'{I),P''{fr) + P''{I)] = [P'',P' 

= 0, 
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which means the total four-momentum operator, the sum of free and inter- 
acting four-momentum operators, also satisfies the point form equations. 

As stated in the introduction, the trilinear vertex is assumed to be bilin- 
ear in fermion-antifermion creation and annihilation operators, and linear in 
boson creation and annihilation operators. That is, such vertices have the 
general form V (a^ + b) (a + b^) (c+c') ^ (a^a + bb^ + a^b^ + ba) (c + c^) so 
the interacting four-momentum operator has the general form 

dvidv2dv{X^^{k)i^i^a\^ai^Ck + X!^2{k)hi2hA2^k 

+X'l2{k)hi2<^\A^Ck + X^^{k)i^i^bi^ai^Ck + he); 
here Xii{k)i-^i^ = F^{vi — V2 — Kv)Mu{iii2k) and 

F'^{u) := / d'^xS{x ■ X — 7-2)a;Me»^-" comes from locality. Mll(^l^2A;) are spinor 
and internal symmetry matrices coming from the free fields, dependent on 
how the fermions and antifermions are coupled to the bosons. The other 
three X's have a similar form. 

Now all quantities except the coupling constant are dimensionless, so g 
has the dimensions of mass. Write g = ma and divide by m. Then all 
four-momentum operators are dimensionless. 

Further, all fermion-antifermion terms are bilinears with the common 
factor Ck or cj[.. So write 

P'{I) = aY.j dv{A{X^,)cu + A{X^,)^cl), (13) 

where A{X'^) := {a\^Mi){X'^)hi^{ai2A2Y- 

In these variables, the eigenvalue problem, Eq.(3), to be solved for trilinear 
interactions is 

{PpiM + Y. j dv{Kv^cick 
+aA{X^,)c, + aA{X^,)^cimx> = A'^I^a > (14) 

This eigenvalue equation is covariant, so Lorentz covariance can be used 
to write = (A, 0,0,0). Then the zeroth component of the eigenvector 
equation is 

+aA{Xl)ck + aA{Xl)^cl))\^x> = A|*a > (15) 
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In particular the vacuum eigenvalue problem is characterized by P°|f2 >= 
and U\\Q >= \Q >. If the point form Eq.(2) holds, then 

UaP^I^ > = UAP^Uj^^UAin > 
= 0, 

which implies that the momentum operator acting on the physical vacuum 
also gives zero, as required. Thus it suffices to analyze the eigenvalue prob- 
lem, Eq.(3) only for the zero component, /i = 0. 

A fundamental difficulty in solving the eigenvalue problem, Eq.(15), is 
that the interaction term takes elements out of the Fock space. To remedy 
this problem, orthonormal bases in the one particle spaces are truncated so 
there are N basic modes, which include the finite velocity modes, as well 
as spin projection and internal symmetry modes, and as N gets larger, the 
continuum (inductive [6] ) limit is approached. 

To distinguish between the continuum energy operator in Eq.(15) and its 
finite approximation, the fundamental operator for trilinear interactions will 
be (a Hamiltonian) denoted by H, made out of creation and annihilation 
operators with a finite number of modes, whose form mimics Eq.(15): 

H = Y.('Mi(^i + blbi + i^clck) + aY,AXk)ck + A{Xl)cl (16) 

= E + E - bibl + kcJq) + a{A{Xi)ci + A{Xj)cl) (17) 

= Ee^ + ^(^) + '^Ee4c^ + aE('^Wc. + ^(^/)4), (18) 
E : = (iza5f(ei, 62,..., Cat, -61,-62,..., -Cat), (19) 

where the discrete "energy" = \Jl + vf, and the A{) notation emphasizes 
the unitary structure developed in the next section. 

3 U(2N) Fermionic Structure 

The fermion-antifermion creation and annihilation operators appear in the 
Hamiltonian, Eq.(18), only in the combinations a|aj, 6^6], a|6] and bittj. In- 
dices range from 1 to N (number of fermion modes). Thus, consider the 
correspondence a] ^ A}, Ui ^ Ai,bi ^ A+at, where the indices 
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on the A's range between 1 and 2N. Then the bihnears A^A generate the Lie 
algebra of U{2N), with commutation relations 



[AlAf3, AlA,] = AiAJp, 




(20) 



For example 




A\j^^Aij^N5j^k - A\A.j5i^i 



as required. 

Again use the definition A{X) := A^^X^pAp,! < a, f3 < 2N; then all 
fermionic terms in H are of this form; for example, the free fermionic energy 
is 



It should be noted that bb^ and not b^b is the proper element in the Lie 
algebra of U{2N), Eq.(20). Further, there is no longer any trace of fermion 
and antifermion anticommutation relations; all of the fermionic structure is 
given by the Lie algebra of U{2N). Only the representations of this Lie 
algebra can be traced back to the underlying fermionic structure. 

The representations of U{2N), in Gelfand-Zetlin labeling [9J, are all the 
antisymmetric representations, written (!,...,!, 0,...0), of length 2N, with the 
identity representations given by all zeroes or all ones. The antisymmetric 
Fock space is the direct sum of all these irreducible representation spaces, 
and is of dimension 2^^. 

But a more convenient way of labelling the fermionic representation spaces 
is with the "baryon number" operator. Define B + N := + ^i^l) = 

J2 -^{1)^ a first order Casimir operator with eigenvalues from -N to +N. Each 
integer corresponds to a given irreducible representation. For example the 
B = sector corresponds to the Gelfand-Zetlin label with equal numbers of 
zeroes and ones. 

In the B = sector, fermionic states can be written as \^i^2---^^n >, where 
the //'s form a shuffle ranging between 1 and 2N[T0]; a concrete realization 
is given by minors of determinants, A^'^'"^(z), with z a 2Nx 2N matrix[n]. 



Hpifr) 



kbl) 
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Writing states as |/ii/i2.../UAr > is equivalent to writing al^bj^...al^bj JO >; in 
this notation the Fock vacuum is |0 >= + 1...2N >~ A]x^':j:'];;;^(2;), while 
the highest state is al^b]^...al^b]jO >= \1...N >~ A\---^{z). 

Since fermion-antifermion operators appear only as bilinears in the Hamil- 
tonian, Eq.(18), a fundamental operator identity which makes such a basis 
representation useful is 

AiAi3\fii...fiN > = 5^1/3 1 a... /i AT > +...6f,^Jfii...a > (21) 

All other baryon number subspaces can be written as |/ii/i2---/ifc >,k = 
1...2N, and an operator identity similar to Eq.(21) applies also to these sub- 
spaces. 

To get from the B=0 subspace to other baryon number subspaces, it 
is possible to use products of operators, A'^^_^...Aj^^. Though single operators 
anticommute among themselves, they satisfy commutation relations with the 
bilinears: 

[AiA^,4] = dp.Ai (22) 
[AiY),Al] = Y^,Ai (23) 
[AiAp.AJ = -5^^A0. (24) 

Here the contrast can be seen with the representation structure of the 
usual fermion anticommutation relations. The fermionic irreducible repre- 
sentation spaces of U(2N) are given by the baryon number operator eigen- 
values, ranging between -N and +N. Since the Hamiltonian is bilinear in 
fermion-antifermion operators, it never mixes these spaces. However, single 
fermion operators do mix these spaces, resulting in there being only one ir- 
reducible representation for the fermion commutation relations, which is the 
direct sum of all the baryon number subspaces. In the limit when N goes 
to infinity, the representation structure radically changes; it is this inductive 
limit [B] that will be studied in the next paper. 

4 Bosonic Contractions 

Even when the number of modes N is finite, it is still possible to have an 
infinite number of bosons in each mode. In this section we show how to con- 
tract a unitary Lie algebra [TT] to the full boson algebra; since the irreducible 
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representations of a unitary algebra are finite dimensional, the contraction 
is done in such a way that as the contraction parameter goes to zero, the 
irreducible representation goes to infinity. 

Consider then the bosonic Lie algebra consisting of creation and annihi- 
lation operators, c|,Cj. Adjoin to this the elements Lij := clcj, so that the 
commutation relations of the four elements are 

[L^J,ci] = cl6,k (25) 

[Lij,Ck] = -CiSjk (26) 

[hjj] = (27) 

h,c]] = Ui, (28) 

[q,/] = (29) 

The commutation relations are those of the semidirect product of the unitary 
algebra with the Heisenberg algebra. 
For each mode this gives 

[L,c^] = (30) 
[L,c] = -c (31) 
[c,ct] = /, (32) 

which is the usual harmonic oscillator algebra for each mode. 

Consider next a U(2) Lie algebra, with elements Ji, J2, and J± and the 
following commutation relations: 

[Ji,Ji] = (33) 

[Ji,J±] = ±J± (34) 

[J2,J±] = {-)±J± (35) 

[J_,J+] = J2-J1 (36) 

Now modify the basis of this Lie algebra by defining J± := pj± and 
J2 := p^J2, with p a positive number; then the Lie algebra becomes 

[Ji,J2] = (37) 
[Ji,J±] = ±J± (38) 
[J2,J±] = ±{-p')J± (39) 
[J_,J+] = J2-PV1 (40) 

(41) 



11 



In the contraction limit when p — this Lie algebra agrees with one mode 
bosonic algebra (J2 plays the role of the identity operator) 

Next consider a concrete realization of the U(2) Lie algebra, with 

J, ^ 3| (42) 

J, - -I, («) 

■U - =^ (44) 

J. - 4 (45) 

on the holomorphic Hilbert space of two complex variables[9]. Then the 
bosonic representations, labelled (M,0), are the homogeneous polynomials of 
degree M, with an orthonormal basis given by < wlM, n >= " = . In 

A/n!(M— n)! 

order that such a basis be holomorphic, the total number of bosons in any 
mode is restricted by M. Further 

< z,w\J+\M,n> = — ' 



^n\{M -n)\ 



= Vn+l^Jp^{M -n) < z,w\M,n + l> . (46) 

When M — > cxD, p ^ such that Mp^ = 1, the usual boson calculus result is 
recovered. 

This can all be generalized to N modes; now the Lie algebra is U(N+1) 
with generators 

J, - 3.^ (47) 

Jt ^ (48) 
J- ^ .A (49) 

- (50) 

As with one mode, pJj~ — > q, J2 — > as p — * 0. 
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Bosonic representations of U(N+1) are written as (M,0,...0)[n]; polyno- 
mial realizations of basis states are given by 



< z, w\Mn > 



(51) 



Jj I M, max > 



Jii\Mn > 
J2\Mn > 

j-\m6> 



0. 



ni\Mn > 







{M -Y,n^)\Mn> 



(52) 
(53) 
(54) 
(55) 



The idea now is to replace a Hamiltonian given in terms of boson cre- 
ation and annihilation operators with the J operators. For example, the 
Hamiltonian, eq.(17), is replaced by 



Hm = Eei + Ee..(ak-Ml + «:^n) + «(^(X,)J-+^(X/)J+)(56) 



and in the contraction limit the eigenvalues of Hm should pass over to the 
eigenvalues of H. An example of an exactly solvable model without fermions 
is given in Appendix A. 



5 Modeling the QCD Vacuum: A Simple One 
Mode Model 



Up to this point only trilinear interactions have been discussed in this paper. 
But, as shown in Appendix B, the ground state for Hamiltonians of the 
form given in Eq.(18) as a function of the bare coupling constant a goes 
as -\constant\a; as discussed in the conclusion, it is not possible to add an 
arbitrary constant to the Hamiltonian, so there is no zero energy ground 
state for nonzero a. This means it is necessary to add other terms to the 
model Hamiltonian. One of the goals of this series of papers is to examine 
inductive limits of the QCD Hamiltonian. Thus it is natural to add boson 
self-coupling terms to the Hamiltonian with trilinear interactions. 

But before examining a model Hamiltonian with both trilinear and boson 
self-coupling terms, it is necessary to investigate the convergence rate of 
boson contractions. The previous section showed that boson contraction 
converges to the full infinite dimensional limit. What is not clear is how fast 
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boson contraction converges. A way to study this question is to examine the 
quantum anharmonic oscillator as a one mode field theory. 
Thus, consider the Hamiltonian 



H 



c^c + aHc + c^Y 



(57) 



The finite number of bosons Hamiltonian is given by generators of the 
U(N+1)=U(2) Lie algebra: 



H 



M 



ji + a\j+ + .Ly 

2 

M+ J, 



2 



4 
z 1 



(58) 
(59) 

(60) 



and in the contraction limit, the eigenvalues of Hm, Eq.(58) should converge 
to the eigenvalues of H, Eq.(57). In Eq.(59) the Lie algebra basis of U(2) 
has been written in a U(l)xSU(2) Lie algebra basis, and in Eq.(60), the 
contraction parameter has been eliminated by writing p = The goal 
is to numerically find the lowest eigenvalues for fixed coupling as M gets 
large. Using an SU(2) Lie algebra automorphism to interchange the x and 
z generators generates a tridiagonal matrix in the basis gives in Eq.(46). 
Reference [12] shows that the true eigenvalue is approached for M about 100. 

Next consider a one mode problem coupling fermions and bosons, with 
no quartic boson selfcoupling. The Hamiltonians are now 



H 

Hm 
E 



1 + (a^a - hh^) + c^c + a{A{X)c + A{X^)c^) 
1 + A{E) + c^c + a{A{X)c + A{X^y 



.t. 



1 + A{E) 

1 
-1 



Ji + a{A{X)J. 

1 1 
1 1 



A{X^)J^ 



.X 



(61) 
(62) 
(63) 

(64) 



For B = the fermion space is two dimensional (|0 >= |2 >,a^6"''|0 >= 
|1 > in the notation of section 3) and the boson space is M+1 dimensional. 
When Hm is diagonalized, the lowest eigenvalue for small values of M have 
been calculated [13]. The behavior of the lowest eigenvalue for any of the 
M values is linear decreasing with respect to the bare coupling constant. 
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confirming the general behavior for a trilinear coupling that the ground state 
always decreases linearly in a for sufficiently large a. 

Finally, if the boson self-coupling term, the anharmonic term in Eq.(58) 
is added to the Hamiltonian, Eq.(63), we have a one mode model of trilinear 
coupling with a quartic boson self-coupling, a simple QCD one mode model: 

H^f^'' = l + AiE) + J, + a{AiX)L+AiX^)J+) 

W{J- + J+f] (65) 

the lowest eigenvalue for several small values of M, as a function of the bare 
coupling parameter have also been computed [13J. For odd values of M, the 
ground state eigenvalue as a function of the bare coupling parameter passes 
through zero (when M is even there are spurious zeroes) ; this raises the 
possibility of fine tuning the bare coupling parameter, as discussed in the 
conclusion. 



6 Fermion-Boson Isospin Model 

In this section we analyze the ground state structure of a few mode system, 
generated by isospin internal symmetry and only one velocity mode. The 
system will consist of isospin 1 "pions" coupled to isospin 1/2 "nucleons" 
and " antinucleons" . The Hamiltonian will have "kinetic energy", trilinear 
coupling, and finally, "pion self-coupling" terms. 

To begin, the bosonic isospin operators are generated from a U(N+1)=U(4) 
bosonic symmetry, enumerated as 1,0,-1: 

/f = J„-J_,.,(=.,|--...^) (66) 



rB 



which then satisfy the isospin algebra, [I^ , If] = If, and [If, I^] = 21 f. 
The bosonic states in the U(4) algebra in a holomorphic basis are given 

nj- no M-y^ni. 

as M, n I , no, n_ >= / =; the no-pion state is given by M=0. 

Y'n+!no!n_!{A/-2^nfe)! 

The M=l representation is four-dimensional, with one no-pion state (isospin 
0) and three one-pion states (isospin 1): 
|M = 1, 000 >= w, isospin 
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\M = 1, 100 >= zi, isospin 1, (+1), 
\M = 1,010 >= zo, isospin 1, (0), 
\M = 1,001 >= isospin 1, (-1). 

The M=2 representation is 10 dimensional, with isospin 2 and 1 states, 
and two isospin states; only the and 1 states are listed here: 
|M = 2,/ = 0>=^, 

|M = 2,J = 0>=^ziz_i-^4 
\M — 2,1—1 >— wzi, wzq, wz^i 



More generally the isospin states can be written as 
|/ = >= J2m (^m\M — 2n, n — rn, 2m, n — m > 

In a similar fashion the proton and neutron give 2 modes, which, along 
with their antiparticles, give the unitary group U(2N)=U(4). In terms of 
creation and annihilation operators the four particle states can be written as 

\p> = a\\0 >,\p >= bl\0 >, 

\n> = al\0 > , \n >= bl\0 > . 

Then the fermion isospin structure has the form: 



= -{a[ai 



+ hb\-b2bl) 



J (4^1 -4^2 + 4^3-^1^4 



2 
A 



T3 


" 


1 


' 1 " 





T3 


'^' = 2 


-1 



AU2 + 4^4 

[A 



A 










' 


1 " 
















T3 





] = [A 









7"3 





,A 



,A 



r+ 


r_ 









27^ 



(68) 
(69) 

(70) 

(71) 

(72) 



and the action of the operators on the primitive states is given by 

P >— ||P >)-^3^l^ >— "11*^ > B=l and by 
>, /fin >= ||n > for B=-l. 



-'3 
-'3 



The B=0 space is 6 dimensional, with elements like |0 >,a|62|0 > .... 



In A^A^ notation, |0 >^ |3,4 >,al6^|0 >^ |3, 1 > 
a[b{albl\0 >^ |1,2 > .... 
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There are 3 isospin states: 



|/ = 0>i = |3,4>(=|0>) 

|/ = 0>2 = \l,2 > {= a\b\albl\0 >) 

l^ = 0>3 = -^(|1,4>-||2,3>) 

(=;^(al6l|0> -46^10 >) 



Given the boson and fermion Hilbert spaces, we wish to construct a 
fermion-boson trihnear couphng Hj,, which is an isospin invariant: 



I3 = I[ + 1^ = A® I + I ® ( Jii - J_i,_i) 
/+ = I^ + I^ = A(^I + I® \/2( Ji,o + Jo,-i 

Hj = aJ2iA{X,)0j^+A{Xl)0j^), 

with the matrices Xk,Xl determined by [h,!!] — [I+,H] — 0; 
this gives 



(73) 
(74) 

(75) 

(76) 



Xi = -Xl,-- 
Xo = Xl^-V2 

X_i = -xl^ 



n T3 

T- T- 
T- T- 



and the interacting Hamiltonian has the form 
Hi 



A{Xi) ( Jf + Jli) + AiXo) ( Jo+ - Jo) 
+AiX^i) ® ( Jli + Ji+) 
(4 + 4)(A2 + ^4)®(Jr- J1i) 
--^(Al + 4)(Ai + A3) ® (Jo- + 4) 



(77) 
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+-^(4 + 4)(^2 + A4) ® (Jo + 4) 

+{Al + 4)(Ai + A3) ® {-JZ, + A). (78) 

Given the isospin Hamiltonian, Eqs.(75),(76), we wish to solve the "vac- 
uum" problem H\^ >= XminV^ > for the "trilinear interaction" on the B=0, 
1=0 space. M=l generates 4 isospin states; in particular 



|J — >4 — ^Cq]^_J\I 



1 



|1,3> 



-1 



l,m>F 
1 

72 



+ |2,4> 
i//(|3,4 > ®w) 
V3|/ = 0>4, 
Hi{\l,2>®w) 
-a/3|7 = 0>4, 

i/,-^(|2,3>-|l,4>)®«;) 

= 0, 

Hi\I^Q>i = V3(|/ = >i -|/ = >2). 
Then the Hamiltonian matrix to be diagonalized is 

-2 x/3a 



1=1, —m >B 
(|l,4>+|2,3>)®^o 



i/7|/ = 0>i 

Hj\I = >,s 



i7 



2 -\/3a 

y/Sa —\/3a 



(79) 



and has lowest eigenvalue Xmin{(^) = — V4 + 60;^. For a ^ 1 the eigenvalue 
goes as a negative constant times a, consistent with the behavior of trilinear 
couplings discussed in Appendix B. A factor 2 should be added to the free 
Hamiltonian, eq.(75); then the minimum eigenvalue is zero when the bare 
coupling is zero. 

To model QCD behavior a quartic boson self-coupling term should be 
added. Though there is not a sufficiently rich isospin quartic operator that 
mimics the anharmonic oscillator, a possibility is 



H, 



self 



a 



^J+J,-|/ = 0>4 = 3|J=0>4, 
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so the Hamiltonian to be diagonalized now is 



H 



Ho + Hj + Hself 

-2 
2 




(80) 











(81) 



in this simple model the "self-coupling" is not sufficient to produce an eigen- 
value zero for the lowest eigenvalue. In the next paper models with many 
velocity modes will be analyzed, adjoined with an internal symmetry of the 
type analyzed in this section, to see if "self-coupling" is able to produce a 
zero eigenvalue for a particular value of the bare coupling constant. 



This paper is the ffist in a series of papers that explore the possibility of ap- 
proximating solutions of quantum field theories through the use of inductive 
[6| and contractive limits [11]. The starting point is an algebra of operators 
bilinear in fermion and antifermion creation and annihilation operators which 
generate a unitary Lie algebra. For infinite degree of freedom systems there 
is a very large class of representations of this unitary algebra. The idea is to 
explore the inductive limit of nested finite dimensional unitary subalgebras 
of the full infinite dimensional algebra. 

But before this inductive limit can be explored it is ffist necessary to deal 
with the boson algebra. Here also there is an algebra of operators which, for 
infinite degree of freedom systems, have a rich representation structure that 
can be explored as an inductive limit. The problem that aries for bosons is 
that, even before considering an inductive limit, it is necessary to deal with 
the fact that in any mode, there can be indefinitely large numbers of bosons. 
One of the main goals of this paper has been to show how contraction limits 
provide a natural way of approximating indefinitely large numbers of bosons. 
That is, if the bosonic Lie algebra, the semi-direct product of the Heisenberg 
group with a unitary group (see Eqs.(25) through (29)), is replaced by the 
compact group U(N-|-1), all of the irreps, and in particular the bosonic irreps 
of the form (M, 0,...,0)[9], are finite dimensional. Here M is an irrep label that 
specifies the maximum number of allowed bosons in all of the modes; as seen 



7 Conclusion 
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in Eq.(55), the raising operator, annihilates the state with the maximum 
number of bosons. This is to be contrasted with simply putting a cutoff on 
the maximum number of allowed bosons, for here it is the Lie algebra itself 
that keeps the number of bosons finite. When the contraction parameter p 
goes to zero as the irrep label M goes to infinity, such that Mp^=l, as shown 
in section 4, the full boson algebra is recovered. Appendix A shows how the 
contractive limit is obtained for a simple exactly solvable model. Further, a 
simple one-mode model based on the quantum anharmonic oscillator shows 
that full convergence to the limit for the ground state already occurs for M 
about 100; more work is required to see how large M must be for systems of 
many modes. 

The philosophy underlying these papers is to preserve as many symme- 
tries as possible, even when the underlying variables are discrete and range 
over finite values. Symmetries here include not only the usual Lorentz and 
internal symmetries, but also the unitary symmetries leading to inductive 
limits generated by fermion-antifermion bilinears (in creation and annihila- 
tion operators) and boson symmetries, with both inductive and contractive 
limits; then interactions are given by products of generators of these Lie alge- 
bras, which in turn generate the four-momentum operator as the fundamental 
operator. 

The context of this paper is point form quantum field theory [3], and the 
basic equations to be satisfied are the point form equations (1) and (2). 
In the point form all interactions are in the four-momentum operator and 
the Lorentz generators are kinematic. Eigenvalues of the four-momentum 
operator, Eq. (3) generate the observables, the physical vacuum, the mass 
spectrum, and scattering states. Lorentz transformations are automorphisms 
on the algebra of fermion and boson operators; when the algebra is finite di- 
mensional (for a finite number of velocity modes) the Lorentz transformations 
are realized as a finite subgroup of the permutation group. 

A second goal of this paper has been to study generic properties of Hamil- 
tonians with trilinear coupling, such as Eq.(16). The motivation here is to 
study the vacuum properties of such Hamiltonians, in preparation for a study 
of the vacuum structure of the full four-momentum operator. A vacuum so- 
lution to the eigenvalue equation (3) means that in a suitable representation 
space of the algebra of operators, the vector \Vt > must satisfy 







n >; 



(82) 
(83) 
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what is of interest here is that if the point form equations (1) and (2) are 
satisfied, it is not possible to add constants to the four-momentum operator 
and still satisfy Lorentz covariance, Eq.(2). Now the free four- momentum op- 
erator satisfies the above vacuum equations, where the vacuum is the usual 
Fock vacuum. As shown in Appendix B, the generic structure of Hamiltoni- 
ans with trilinear couplings is that the vacuum energy decreases linearly with 
the bare coupling constant. This means there can be no solution to the vac- 
uum problem, Eqs.(82),(83), with just trilinear couplings. Other interactions 
such as a boson self-interaction are necessary for the vacuum solution to give 
zero when acted on by the four-momentum operator. Section 5 presented a 
simple one-mode model in which the boson self-interaction is given by the 
quantum anharmonic oscillator. When combined with a trilinear coupling, 
the eigenvalue problem for small values of M gives nontrivial solutions to the 
vacuum problem. This raises the possibility of fine tuning the bare coupling 
constant, though, of course, much work remains to see whether this sort of 
fine tuning persists in the inductive limit. 

The longer term goal of this series of papers is to investigate the sorts 
of problems raised above when the number of velocity modes gets large, 
and in particular to investigate the nature of the physical vacuum that is 
generated by the QCD Hamiltonian. Here also there will be an interplay 
between trilinear interactions, of the quarks with the gluons, and gluon self- 
interaction terms. 

Acknowledgement: I wish to thank Fritz Coester for many stimulating 
conversations. I also thank Addison Stark and Kevin Murphy for results and 
discussions on the anharmonic oscillator and one-mode trilinear coupling 
Hamiltonians. 



8 Appendix A: Exactly Solvable Bosonic Con- 
traction Model 

Consider the exactly solvable boson Hamiltonian, 

H = ^e4Q + a^(AQ + A*4) (84) 

= E^^(4 + ^)(Q + ^)-«'E^; (85) 

\{gnd) = -a^^^. (86) 
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In Eq.(84) the fermion operators in Eq.(18) are replaced by constants Dj. As 
is easily checked, Cj + also satisfies bosonic commutation relations with 

its adjoint, and hence, the spectrum of is ejnj — a^^^^; the ground state 
is given by all the n's zero, which is Eq.(86). 

The Hamiltonian Hm which contracts to the Hamiltonian H of Eq.(84) 
is made out of U(N+1) Lie algebra elements; the goal of this appendix is to 
show that the ground state of Hm contracts to the ground state of H: 



Hm = Ee^^- + "E(A^r + A^^/i") (87) 

= TrMJ (88) 

= TrU^AUJ (89) 

= TrAUJU^ (90) 

= ^A,j; + A^+ij;; (91) 

here the matrix M = ( ^'^Z ' 1 and J = f =U^j'U 

is an automorphism of the U(N+1) Lie algebra. 

Corresponding to the change of the Lie algebra basis is a change of the 
representation basis: 

<z,w\Mn>' = < {z,tv)U\Mn>; (92) 
HM\Mn>' = (5^ A,n, + A^+i (M-^n,)) I Mn>' (93) 

= (^(Ai-A7v+i)ni + Ajv+iM)|Mn>'; (94) 

Hm\MQ>' = MAjv+i|MO>' (95) 

is the ground state ( with eigenvalues ordered by Ai > ... > Ajv > Ajv+i). 
It is simplest to compute the contraction limit when there is only one 

mode (N=l); then = ^ ^ with eigenvalues A± = ^^^^ +4a _\D\_^ 

The contraction limit is obtained by replacing a with ap and letting p — > 
0, M — > 00 such that Mp^ = 1. Inserting the contraction parameter into 

the miniumum eigenvalue, Eq.(95), gives MA_ = ^(1 - y'TT^^^V) ^ 

— " as the ground state eigenvalue in the infinite boson limit. This result 
agrees with the ground state, Eq.(86) for one mode. 
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9 Appendix B: The Ground State Eigenvalue 
in the Strong Couphng Limit 

The goal of this appendix is to show that the Hamiltonian for trihnear 
fermion-boson interactions has the property that the ground state goes as 
a hnear function of the bare couphng constant a, for a ^ 1. Moreover, the 
constant multiplying a is negative. The general form of the Hamiltonian 
is given in Eq.(18); here it is assumed that the matrices commute for 
different modes, [X^,X^] = 0: 

e,(ala, + bib, + kc\c,) +aY,A{X^,)c, + A{X^,)^c\- (96) 

+- E + + K,c\ci) (97) 

a 

E(4 + ^(X°))(q + ^(X°)t) - Y.Ax^)Ax^)' (98) 

- E + - E (ctlcii + b\bi + hic\ci). (99) 
a 

The part of the Hamiltonian given in Eq.(98) has a ground state eigenvector 
and eigenvalue given by 

1*0 > = 7Ve-S'^(^?)^^l|A> (100) 
^(X°)t|*o> = -Afe|*o>, (101) 

where is a normalization factor such that < \E'o|*o >= 1- Further, 
Cfel^^o >= -A^l^o >, so that < ^o\A{Xl)A{Xiy + 4cfc|*o 2|Afe|2. 
Finally, the part of the Hamiltonian multiplying - in Eq.(97) is a nonneg- 

ative operator, < \I^o| E ei(aja, + b\b, + kc|q)|*o >:=< *o|^o|^o >> 0, 
so that, in first order perturbation theory, the lowest eigenvalue of is 
— 2X] |Ap + <'^o\Ho\^(,> ■ then follows that the lowest eigenvalue of H in 
first order perturbation theory is — 2q:X1|AP+ < \I/o|i^o|^o >• That is, as 
a function of the bare couphng constant a, it goes linearly with a, with a 
negative slope. 



H = 

H 

a 
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